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We consider the one-dimensional Glauber dynamics with coupling disorder in terms of bilinear 
fermion Hamiltonians. Dynamic exponents embodied in the spectrum gap of these latter are evalu- 
ated numerically by averaging over both binary and Gaussian disorder realizations. In the first case, 
these exponents are found to follow the non-universal values of those of plain dimerized chains. In 
the second situation their values are still non-universal and sub-diffusive below a critical variance 
above which, however, the relaxation time is suggested to grow as a stretched exponential of the 
equilibrium correlation length. 

PACS numbers: 02.50.-r, 05.50.+q, 75.78.Fg, 64.60.Ht 

I. INTRODUCTION 

In the absence of general principles accounting for the evolution of non-equilibrium systems, kinetic Ising models 
have been used to characterize a large variety of slow non-equilibrium processes [l], [|| ■ Among the most ubiquitous 
of these latter, the coarsening dynamics after a sudden quench from a disordered phase at high temperatures to an 
ordered one below the critical point still remains a subject of much interest and debate For pure substrates it 
is by now established that coarsening domains are characterized by a single time-dependent length scale growing as 
t 1 / 2 , where z is the dynamic exponent of the universality class to which the dynamics belongs. Depending on whether 
or not the order parameter is preserved by this latter, e.g. Kawasaki Q or Glauber Q dynamics, the network of 
(ferromagnetic) domains will coarsen respectively according to a sub-diffusive Lifshitz-Slyozov growth (z = 3, see 
however Ref.p ), or following a plain diffusive Allen-Cahn behavior [z = 2) 

On the other hand, the case of coarsening systems with quenched disorder substrates modelling the effects of various 
experimental situations has also received considerable attention [H, 0] , especially in the context of random field and 
random bond Ising models [8l4l2|, In general, it is expected that quenched impurities play the role of energy barriers 
slowing down coarsening domains, the boundaries of which move by thermal activation over the landscape of these 
trapping centers 0]. However, it is not yet clear whether universal growth can take place at late evolution stages in the 
presence of strong disorder. Recent Monte Carlo simulations [l2| of the Glauber dynamics using random ferromagnetic 
couplings drawn from uniform probability distributions suggest that after long pre-asymptotic crossovers the final 
growth in one dimension (ID) recovers the usual diffusive behavior referred to above, while becoming logarithmically 
slow in the two-dimensional case. For the ID situation this is somehow intriguing as already at the level of a plain 
alternating-bond or dimerized chain [l3| the dynamic exponents are known to be non-universal and sub-diffusive 

As part of the ongoing efforts in this context, here we revisit the ID Glauber disordered dynamics in terms of bilinear 
fermion fields associated to the kinetics of its kinks or domain walls. Here, we rather focus on the dynamic exponent 
characterizing the actual growth of relaxation times r with equilibrium correlation lengths £, which according to 
critical dynamic theories [Ty| should scale as r cx £ x . But it is known 0, [ljj that this latter exponent coincides with 
z because at times comparable with r the average domain size becomes of order £ , so both descriptions are ultimately 
equivalent at large scales of space and time. The idea is to avoid the problem of dealing with prohibitively long 
transient regimes by pinpointing directly the relaxation time embodied in the spectrum gap of the evolution operator 
which, for that purpose, will be constructed and diagonalized in a domain wall representation. Although the critical 
point is strictly zero, the analysis is kept within low temperature regimes (with £ and r becoming both arbitrarily 
large there), otherwise the dynamics would be rapidly arrested owing to the proliferation of large quantities of 
metastable states. The disorder considered throughout may include either exchange couplings of a single type (ferro 
or antifcrromagnetic) or mixed ones, though as we shall see in a moment, the dynamics of these situations can be 
mapped onto each other via a simple spin transformation. 

The outline of this work is organized as follows. In Sec. II we describe the basic kinetic steps under random spin 
exchanges using a kink or dual representation, and write down the underlying Glauber generator in quantum spin 
notation. Exploiting detailed balance [l5j we then recast this latter operator in terms of bilinear and symmetric 
fermion forms, thus allowing to evaluate its spectra through the solution of a secular problem whose dimensions grow 
linearly with the system size. In Sec. Ill we present the results arising from the numerical diagonalizations of this 
problem both for binary and Gaussian disorder realizations. Non-universal dynamic exponents as well as scaling 
forms drawn from the exact solution of the dimerized chain are then proposed and seen to be consistent with some of 
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our numerical findings. However, when the width of the Gaussian disorder exceeds a certain threshold the relaxation 
time grows in a rather different (non-power-law) and much faster form. We close with Sec. IV which contains a 
recapitulation along with brief remarks on open issues and possible extensions of this work. 



II. DYNAMICS AND BILINEAR FORMS 



Let us consider the Glauber dynamics of the Ising chain Hamiltonian % — — ^2n Ji Si Si+i with L spins Si = ±1 and 
disordered couplings Ji chosen arbitrarily with any sign and strength. Unless stated otherwise, periodic boundary 
conditions (PBC) will be used throughout. The system is in contact with a heat bath at temperature T, which 
causes the states = \S\, ... ,Sl) to change by random flipping of single spins. The probability per unit time 
W(Si — > —Si) or transition rate for the i-th. spin to flip is set to satisfy detailed balance [l5[ (see below), so that the 
system relaxes to equilibrium. As is known, the simplest choice of rates complying with this condition corresponds to 
that of Glauber la 



1 - f (% Si-i +7+^ +1 ) 



W(S^-S % ) = 

7 ± = tanh {Ki + K^) ± tanh [K t - , (1) 



where Ki = Ji/ksT, and the spin flip rate a is hereafter taken as 1/2. We shall always measure the temperature in 
energy units; equivalently the Boltzmann constant ks is set to 1. Note that the nature of the dynamics is basically 
unaltered by whether some or all couplings are ferromagnetic (Ji > 0) or antiferromagnetic (Ji < 0). To check out 
rapidly this issue consider for simplicity a chain with open boundary conditions (just to avoid frustration due to 
eventual mismatch of antiferro exchanges), along with the mapping 

i-l 

Ri = Si , R, = Si J] sign (J n ) for 1< i < L , (2) 

n=l 

to new spins Ri = ±1. Then, to preserve the energy of the mapped configurations and therefore to maintain invariant 
all transition rates (these in turn depending on AH /T) , clearly in the equivalent i?-system all couplings must become 
ferromagnetic because Si S^+i — > sign ( Ji) Ri Ri+\. In particular, the dynamics of the random ± J chain would then 
reduce to the homogeneous one. Also, notice that the single-spin flip dynamics maps onto itself. In that latter respect, 
this argumentation would not apply to the Kawasaki dynamics [ij as owing to the exchange of S-pairs under antiferro 
couplings, the mapping ([2]) would then allow parallel i?-pairs to flip. 

Strictly, Ising Hamiltonians possess no intrinsic dynamics since all involved spin operators commute with one 
another. But when these systems are endowed with extrinsic transition rates W(S — > S"), such as those of Eq. (TTJ), 
their time evolution is described in terms of a Markovian process governed by a master equation [151 ] . For our 
subsequent discussion it is convenient to think of this latter as a Schrodinger equation in an imaginary time 

d t \P(t)) = -H\P(t)), (3) 

under a pseudo Hamiltonian or evolution operator H. This provides the probability distribution to find the system 
in a state \P(t)) = J^s ^(^i 1^) a ^ time t, from the action of H on a given initial distribution, i.e. \P(t) ) = 
e -Ht |p(o) ). To ensure equilibrium at large times, the detailed balance condition referred to above simply requires 
that P eq (S) W(S S') = Peq(S') W(S' S), V \S), \S'). As usual, the diagonal and non-diagonal matrix elements 
of this Markovian operator are given respectively by 

(S\H d \S)=J2 W ( S ^ S ')> {S'\H nd \S) = -W(S->S'). (4) 

S'^S 

In particular, the first non-zero eigenvalue E\ of such stochastic matrix singles out the relaxation time r in which 
we are interested, i.e. r = 1/Rc E\ > 0, which is the largest characteristic time for any observable at a given low 
temperature. The case of Eq — just corresponds to the stationary mode. 

In order to build up and diagonalize the operational analog of Eq. (|4]), in what follows it is convenient to work 
instead with the associated kink or domain wall dynamics rather than with the spin flipping process itself. Both 
representations are depicted schematically in Fig.[TJ Clearly from Eq. (fTJ), the kink pairing and hopping rates indicated 
in this diagram become respectively 

e ±Pi = (l±tanhPi)/2 , P i =K i +K i _ 1 , 

e ±0 . = (l±tanhQi)/2 , Q i = K i -K i - X . (5) 
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FIG. 1: Disordered transition rates [Eq ([5]) ] for (a) creation-annihilation of kink pairs , and (b) hopping of kinks. These latter 
are denoted by vertical lines separating domains of opposite spin orientations. 

Thus, if we think of the corresponding kink occupation numbers (0,1) as being related to usual spin-^ raising and 
lowering operators a + , er~ , we can directly identify the non-diagonal operators H n d = H pa i rs + Hh ops associated to 
Eq. ((U) . Evidently, the transitions of these dual processes will be provided by 

Hpatrs = - ( e -p, a t-i a t + e Pl a l a i-x ) > ( 6 ) 

i 

H hops = - ( € -Qi a i + € Qi a i-l a i ) • ( 7 ) 

i 

On the other hand, a significant simplification arises for the diagonal terms needed for conservation of probability. The 
former basically count the number of hopping and pairing instances at which a given kink configuration can evolve to 
different ones in a single step. Although this would introduce two-body interaction terms of the form &i~ a ^ (7 i+i a i± l 
it can be readily verified that their coefficients will all cancel out so long as the rates involved are those of Eq. ([5| , 
ultimately stemming from detailed balance. Thus, after some brief calculations we obtain 

H d = J2 ( e -n + h '^i a ^) ' 

i 

hi = ( tanhQj + tanh Pi + tanhPj-|_i — tanhQi-|_ij/2 , (8) 

and no kink interactions, irrespective of disorder in the original couplings. 

To recast the non-diagonal parts of Eqs. ((6|) and into a symmetric representation, we have recourse once more 
to detailed balance and rotate each cr^ around the z direction using imaginary angles 4>j — iKj. So, we consider the 
diagonal non-unitary similarity transformation S = exp ( | -fQof ) for which it is straightforward to show that 



of °U ^ at °tx , 



(9) 



besides keeping unaltered the algebra of <r + , a . After this pseudo spin rotation the above non-diagonal operators 
therefore transform as 



Hpairs 



H, 



hops 



~ J2 sechP (<r+ cr+ 1 + H. c. ) , (10) 

i 

-\ sechQ, ( a+ oT_ x + H. c. ) , (11) 



while leaving Eq. © invariant. Thus, we are left with a bilinear Hermitian form whose spectrum gap embodies the 
wanted dynamic exponents referred to in Sec. I. In passing, it is worth mentioning that Eqs. ([5]l. (|10[) and (|lip also 
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define an anisotropic XY chain under both inhomogeneous couplings Ji Y = — (sechQi ± sechP^) and transverse fields 
hi. However, the commutation algebra of its constituents Pauli matrices, as well as that of the above raising and 
lowering operators, would complicate the analysis. Fortunately, since the parity of kinks is conserved throughout and 
all couplings extend just to nearest-neighbors, a Jordan- Wigner transformation to spinless fermions Ci [1 81 ] enables 
us to progress further. With the aid of the real symmetric and antisymmetric tridiagonal L x L matrices (with 
boundaries) A and B, defined respectively as 



A = 



( hi x 2 

x 2 h 2 x 3 
x 3 



V -xi 







-xi \ 



Xl-i 

XL-l h L -i X L 

x L h L 



B = 



( y 2 
-2/2 ?/3 

-2/3 '•• 



V -2/i 







Vi \ 




VL-1 

Vl-i y L 
-Vl / 



(12) 



Xi = — - sechQi 



.!/i 



sech Pi , 



the Jordan- Wigner representation of spins finally lead us to the bilinear fermionic form 



H 



E 



C i A-i,i Cj 



c\ Bi, c 



H.c. 



(13) 



(14) 



The boundary matrix elements in A and B are especially aimed for PBC by respectively taking account of the 



contributions of a\ a 1 — > —c L c\ and er^cr^ 



-> -c L c\. 



Since the number of kinks J2i a i cr i = J2i c l c » under PBC is 



always even, only anticyclic conditions must be used in those boundary terms. 

Following the standard literature [l8j . this quadratic form can be recasted by means of a canonical (unitary) 
transformation into new 6 n -fermions such that [ b n , H 1 = E n b n , and therefore 



H = 



L 



bn - 



(15) 



Here, the additive constants are obtained by the invariance of the trace of H, whereas the eigenvalues E n of its 
elementary excitations are given by the solutions of either of the following two L x L secular equations [l8[ 

(A±B)(A T B)<p± = El<pi. (16) 

Because (A + B) T = A — B the above products are symmetric and hence diagonalizable, both being in turn represented 
by real five-diagonal matrices (plus boundaries), as if they were single particle Hamiltonians with hoppings terms up to 
next-nearest neighbors. On the other hand, since by construction H is a stochastic operator its vacuum energy always 
vanishes and therefore all eigenvalues should be constrained as ^ n E n = L, with E n > 0. This is an issue which we 
shall make use of as a consistency test in the numerical diagonalizations of Sec. III. The wanted relaxation time will 
then emerge by averaging the spectrum gap of this latter secular problem over several of its disorder realizations. 



A. Dimerized case 



Before continuing with the numerical analysis of Eq. (| 16[) we pause to consider briefly the dynamics arising from 
a periodic array of J\ - J 2 couplings over odd and even bonds, say. Equivalently, the above P's and Q's become 
P = K\ + K 2 , Q = K 2 — K\ . As already this situation is capable of yielding non- universal exponents [13| , it is 
instructive to see how these are recovered within the context discussed so far. We begin by Fourier transforming the 
Ci operators of Eq. (fT4")l to a set of wave fermions f q , g q 

L/2 

e^Y^e-^oa, (17) 
4=1 
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with L/2 wave numbers 



qe Q ± = {±2^,±6^,--- ,±(L-6)~,±(L-2)£} , 



(18) 



accounting of the anticyclic conditions referred to above. The e ±lq ^ 4 phase factors are just aimed at producing a real 
representation of H when expressed in terms of /'s and g's. After introducing the 2x2 matrices, 



A„ 



B„ 



parametrized as 



— sech Q cos — , y q 
tanh P ± tanh Q , 



o y q 
y q o 



sech P sin 



2 " 



(19) 



(20) 



it is readily found that the problem splits into subspaces involving only the occupation numbers and pairing states of 
four fermions. Specifically, the dynamics is decomposed as 



H = Le_ p + H A (q) + \ 



(21) 



where 



H A {q) 



q £ 2=" 



qeQ+ 



(ti \ 



B q 
-B q 



9t 

\gl q J 



H.c. 



(22) 



which is the bilinear counterpart of Eq. (| 14|) in this momentum space. Hence, the solution of the secular equation (fTB"]) 
here reduces to the diagonalization of either of the products of the 4x4 matrices 



D„ 



( A q ±B q 
\±B_ q A_ q 

ng the 

ultimately lead to two elementary excitation bands, namely 



(23) 



Thus after straightforward algebraic steps and taking the positive roots of the associated E q eigenvalues, we are 



T ± = 



1 ± y/T+ + T-cosq , 
( tanh 2 P± tanh 2 Q) /2 . 



(24) 



Parity conservation of original kinks allows to create only an even number of these excitations. Then, by occupying 
two E~. levels with q* = ±2n/L when \P\ > \Q\, or alternatively, using q* = ±(L — 2)n/L when \P\ < \Q\, the 
spectrum gap g — r _1 in the thermodynamic limit comes out to be 



.9 = 2(1- tanh-R) , R = max( \P\ , \Q\ ) 



(25) 



On the other hand since in equilibrium all kinks become independent, the spin-spin correlation length is simply 
= — | In ( tanh \Ki\ tanhj^l ) [3, thus growing as £ ~ e 2 '^ 1 ' in the low temperature regime, say for \J%\ > |Ji|. 
Therefore using Eq. (|25]l within that latter limit (g ~ 4 e -2 ^), we can now relate these two seemingly independent 
quantities via the non- universal dynamic exponents of Refs. [131 ]; that is 



z = 1 + |J 2 M| , for |J 2 | > |Ji|. 



(26) 



As expected from the simple but more general arguments given below Eq. @, the role of the coupling signs is irrelevant. 
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Scaling regimes - Finally, and in preparation for the finite-size scaling comparisons of Sec. Ill under both discrete 
and continuous disorder, it is helpful to consider also the spectrum gap of finite chains. Expanding the lower band 
of Eq. (|24|) to second order around either of the q* minima referred to above, and considering afterwards the scaling 
regime L — > oo, £ — > oo while holding x = £/L finite, it can be easily verified that the gap decay exhibits the crossover 

f £~ z for x < 1/tt , 

Y^ 2 ~ z /L 2 otherwise. ^ ' 

Evidently, this implies a typical Arrhenius dependence [2J for the relaxation time at low T's, that is r cx exp(2 6/T') 
for r / 1 ,/i I <C 1, though now with energy barriers b fixed by the interplay between temperatures and sizes, namely 

|J 2 | + |Ji| for T/\J X \ > 1/lnVZ, 

(28) 

I J2 1 — I J\ I otherwise . 

Moreover, both finite-size and finite correlation length corrections to these crossovers can also be estimated by keeping 
next leading order terms in Eq. (|24l) . So, these finite effects are readily found to scale as 

9? = 4 + 4tt 2 (i-^)z 2 ■ (29) 

We will get back to these issues later on in the context of Figs.[3]and|4]of Sec. III. 



III. NUMERICAL RESULTS 



Turning to the diagonalization of Eq. (ITC1) , in what follows we consider the cases of binary and Gaussian probability 
distributions of disorder realizations, both taken site-independent. The coupling concentration {p) of the former is 
such that P(J) = (1 — p) 5j t j t + p6j,j 2 , whereas the latter is parametrized as usual by an average (J) and variance 
a 2 = (J 2 ) — (J) 2 . As the temperature is varied we focus on samples with correlation lengths close to their averaged 
disorder values, so as to minimize the dispersion of relaxation times resulting from finite-size effects. In the discrete 
case this is just equivalent to using a constant number of coupling types, i.e. the correlation length is fixed as 
£ p = — [(1— p) In (tanh \Ki\) +p In (tanh \K^\) ] . For the continuous case, in turn we select and diagonalize only 
those samples having = — ^ m ( tanh \Ki\ ) /L~ — In ( tanh \K\ ), which amounts to work with disorder fixed (or 
nearly fixed) end-to-end spin-spin correlations. 

In Fig. [2] we show the typical probability distributions of gaps arising from this simple procedure. Already for 500 
spins it is seen that standard gap deviations turn out to be fairly small. In particular, the abrupt decrease of these 
distributions above g > (g) possibly signals the presence of natural upper bounds. However, for Gaussian disorders 
with (J I (J) > 0.22 (see below) but holding £ approximately constant (that is to say, choosing a temperature so as to 
keep comparable late domain sizes), this negative skewness is smeared out and the distributions widen significantly. 
As we shall see, at the level of average gaps this latter issue will be reflected in terms of a rather different decay of (g) 
with the correlation length. In averaging these gaps we used chains of up to 3000 spins for which the sampling had 
to be substantially reduced (down to 20 realizations). This is due to, in part, an L 3 increase of the diagonalization 
time of standard routines [19j |. as well as to a slowing down of convergence within low temperature or large £ regimes. 
The smallness of the spectrum gap for these regions and chain lengths also precluded us from using Lanczos-based 
algorithms [2(J, as their convergence became even slower. Ultimately, these issues set the practical limit to the idea 
of singling out arbitrarily large relaxation times, as intended to in Sec. I. 



A. Average binary gaps 



Nonetheless, the relative errors of (g) yet remain small and clear trends emerge as £ is varied. Fig. [3] illustrates 
the resulting average decays for several parameter values of the binary distribution, all with | J2I > |</i|- It can be 
observed that within the very same scaling regime £ < L/tt referred to in Eq. (|27p for the plain dimerized chain, 
our results suggest a similar gap decay i.e. (g) cx \ J ^/ Jl \ and irrespective of disorder concentrations [2ll I22I]. 

Departures of data from this regime are equally followed by the exact solution of the alternating-bond chain, the 
finite-size aspects of which are depicted in the inset. Furthermore, as is shown in Fig.0] above £/L ~ 1/tt the 
numerical diagonalizations also reproduce the universal scaling function contained in the thermodynamic limit of 
Eq. (|29p (appart from an overall p dependent amplitude), which therefore evidence the same crossover decay, i.e. 
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FIG. 2: (Color online) Gap distributions over 3 x 10 5 samples of 500 spins after binning data in 70 intervals using (a) binary 
and, (b) Gaussian disorder realizations with fixed correlation length (£ ~ L/2). For comparison, dashed lines display Gaussian 
fittings with same mean and standard deviation (s) of otherwise skewed distributions, (a) illustrates the typical narrow 
dispersion {s/{g} ~ 0.04) observed in the binary case (here, T/\ Ji| = 0.36, J2/J1 = 2, and p — 0.5). In (b) an almost identical 
quantitative behavior is displayed for (J) = 2 with T / (J) — 0.3 and <r/(J) = 0.15 (innermost curve), whereas by contrast the 
case T/(J) = 0.2, o / (J) — 0.3 > a c (see text) but similar £, develops a more symmetric and broader dispersion (s/(g) ~ 0.12). 
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FIG. 3: (Color online) Average gaps over 20 realizations of 3000 couplings drawn from binary distributions (1 —p) 5j,j x +p5j,j 2 ■ 
Squares, circles, and triangles refer respectively to concentrations p — 1/4, 1/2, and 3/4 with correlation lengths £ p as referred 
to in the text. Solid lines stand from top to bottom for dimerized cases with IJ2/J1I = 1.5, 2, 2.5, 3, and 3.75 (£ = £ p =i/2)- 
Within the scaling region £ < L/tv, these closely follow the non-universal power law decays of their disordered counterparts, 
the gaps of which, for displaying purposes, have been normalized by different amplitudes. The inset illustrates the finite-size 
crossover mentioned in Eq. (|27|l for the dimerized situation using L = 10 3 , 10 4 , 10 5 (top to bottom), and | J2/J1I = 3 (also, see 
Fig-H]). For comparison, the uppermost dashed line in the main panel denotes the uniform case z — 2. 
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FIG. 4: (Color online) Finite-size scaling regimes of average gaps and correlation lengths for a binary coupling disorder with 
p — 0.5 and |Ja/Ji| = 2.5 (main panel). From left to right doted and solid lines stand in turn for L — 3000, 2000, 1000, 500, 
and 200. The data collapse in the scaling region f/L > 0.3 is attained upon setting z = 1 + | JijJ\\. For much smaller values 
of f/L finite-size departures show up with the same sign as those of Eq. (|29|l . the thermodynamic limit of which is denoted by 
the dashed line. Similarly, the inset displays the scaling regimes of a Gaussian distribution with o / (J) = 0.15 upon choosing 
z ~ 2.61. For comparison with the above thermodynamic limit, all finite size data were rescaled by a common factor. 

(g) oc £ 2-z /L 2 , for the binary distribution. Alike the dimerized case, deviations from scaling are negligible within this 
region (all sizes yielding slopes ~ 2), though as £ decreases subdominant l/^ z corrections bring about breakdowns. 
Just as in Eq. (|29[). these happen to be of the same sign and more severe at the smaller sizes. 

B. Average Gaussian gaps 

Similar scaling features appear also for Gaussian disorder provided a is held small enough (see below), although 
finite-size departures become even more prominent for £ < L, as is observed in inset of Fig. 2] In evaluating average 
gaps for this latter type of disorder, and in parallel with the behavior of the gap distributions referred to in Fig.J^b, 
two decay regimes are now obtained according to whether the relative standard deviation a = cr/(J) is smaller or 
greater than a certain threshold a c . This is displayed in Fig.[SJi where we show some trends for 3000 spins using 
several variances for (J) = 1 and 2. As long as a < a c all gaps exhibit usual (but non-universal) power law decays 
oc l/£ 2 with z increasing slightly between 2 and ~2.65 as a increases. The incipient crossovers observed at £ > 1000 
correspond to the scaling regimes already alluded to in the inset of Fig.[4l so there (g) oc £ 2 ~ z . 

By contrast, however, for a > a c a much faster decay emerges over the whole range of accessible correlation lengths, 
whereas the scaling relation put forward in Fig. @] no longer holds. We direct the reader's attention to Fig.[S}D where 
it turns out that the non-linear least-squares fitting of the stretched exponential form 

( 3 )~Aexp(-af), (30) 

is able to follow very closely all numerical diagonalizations. After estimating the stretching exponents of several 
variances it was found that above a threshold of a c ~ 0.22 up to a ~ 0.35, they increase with a as a ~ 6.01 ( a — <r c ) 
as is shown in the inset of Fig.[5p. In nearing that critical variance however, it is difficult to distinguish the power 
law decay obtained before from this new behavior, so the crossover from one cr-regime to the other turns out to be 
smooth. 
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FIG. 5: (Color online) Average spectrum gap for a Gaussian distribution of couplings over 20 samples of 3000 spins, all 
realization having in turn a common correlation length. In the main panels circles and triangles stand respectively for (J) = 1 
and 2. From left to right each data set referrs to (a) a = a/(J) = 0.3, 0.25, 0.2, 0.1 and, (b) a = 0.25, 0.275, 0.2875, 0.3. For 
a < <7 C — 0.22, in (a) data are fitted with non-universal power laws, whereas above such critical value all data closely follow 
stretched exponential decays [Eq. (|30|) ] with stretching exponents a estimated as in the inset of (b). For comparison, dashed 
line in (a) denotes the ordered case a — 0. 



Let us mention that it would be important to elucidate whether this conjectured non-algebraic decay actually 
extends beyond the regions displayed in Fig.[5l>, e.g. £ < 1500, 270 for a = 0.25 and 0.3 respectively. However, above 
those correlation scales the numerical gaps get progressively smaller [ in turn resulting from the even smaller squared 
gaps of the secular problem (116|) ] , which brings about a much slower and erratic convergence of diagonalization 
routines [l9| . 



IV. CONCLUDING DISCUSSION 



Recapitulating, we have constructed a fermionic bilinear representation of the ID Glauber dynamics with nearest- 
neighbor random interactions. This corresponds to the symmetric version [Eq. (J9j) ] of the kink pairing and hopping 
processes schematized in Fig.Q] The relaxation time of these latter was reduced to the evaluation of the single-particle 
spectral gaps of the secular problem given in Eq. (|16l) , the diagonalization of which ultimately led us to non- universal 
forms of domain growth. Some of these were expected and others not foreseen. 

The case of binary disorder lent itself more readily for comparisons with the soluble alternating chain of same 
coupling types (Sec. II A). Although the former situation has no exact solution, based on decimation procedures as 
well as on simple diffusion arguments found in the literature 21] , the numerical coincidence of dynamic exponents with 
those of the soluble case should come as no surprise (Fig. [3]). In addition, the scaling function of Eq. (|29|) reproduced 
our binary data just above regions displaying incipient departures from the actual exponents (Fig. 21 and rightmost 
part of Fig. [3]). Because of this scaling robustness, the same barrier values involved in the crossover of Arrhenius times 
[Eq. (|2"8f ] might be expected to also hold in finite chains with binary disorder regardless of their bond concentrations. 
In passing, it is worth pointing out that such manifestation of finite-size effects has been observed experimentally in 
single-chain magnets [23j. 

Some of these scaling aspects seem to also apply for Gaussian disorder realizations with small relative variances 
(inset of Fig.HJ), for which non- universal but yet power-law forms of gap decay were obtained (uppermost data in 
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Fig. [5^). However, for a j (J) > 0.22 a crossover to new regimes with much faster and non-algebraic decays showed up, 
whilst at the more detailed level of gap distributions this brought about a significant increase of dispersion (Fig.[2}D) . 
The stretched exponential form put forward to fit our data (Fig.Eb) lent further support to the former hypothesis 
which, by inverting Eq. (f5U|) and recalling the large scale equivalences referred to in Sec. I @, [l7|, is tantamount to a 
dynamics of domain scales growing asymptotically as ~ (lnt) 1 /". Curiously, this rather slow form of phase-ordering 
kinetics is reminiscent to that conjectured and studied in late coarsening stages of higher dimensions with bond 
disorder 

Alongside these growth law evaluations, it would be interesting to also address the issue of scaling and superuniver- 
sality (or the lack thereof), in two-time quantities such as autocorrelation and autoresponse functions HHEO]]. That 
research line might well be further investigated with the fermion approach discussed in this work. Contrariwise, the 
effect of external fields 8], whether random or not, would be bound to adopt rather cumbersome expressions had it 
been written in terms of the kink representation. For analogous reasons, the treatment of Glauber dynamics in higher 
dimensions remains beyond the scope of this approach. 
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